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ABSTRACT 

A large-scale study of the molecular clouds toward the Trifid nebula, M20, 
has been made in the J— 2-1 and J— 1-0 transitions of ^^CO and ^^CO. M20 
is ionized predominantly by an 07.5 star HD164492. The study has revealed 
that there are two molecular components at separate velocities peaked toward 
the center of M20 and that their temperatures — 30-50 K as derived by an 
LVG analysis — are significantly higher than the 10 K of their surroundings. 
We identify that the two clouds as the parent clouds of the first generation 
stars in M20. The mass of each cloud is estimated to be 10'^ Mq and their 
separation velocity is ~ 8 km s~^ over ~l-2 pc. We find the total mass of stars 
and molecular gas in M20 is less than ~ 3.2 x 10^ Mq, which is too small by 
an order of magnitude to gravitationally bind the system. We argue that the 
formation of the first generation stars, including the main ionizing 07.5 s tar, 
was triggered by the collision between the two clouds in a short time scale of ~1 
Myrs, a second example alongside Westerlund 2, where a super star cluster may 
have been formed due to cloud-cloud collision triggering. 

Subject headings: ISM: clouds — Radio lines: ISM — open clusters and associations: 
individual: M20 
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Introduction 



High-mass stars play an important role in dynamically agitating and ionizing the 
interstellar medium in galaxies and strongly influence galactic evolution. It is therefore 
crucial to understand their formation mechanisms. Despite numerous theoretical and 
observational studies attempting to address this issue in the last few dec ades, the detailed 



Zinnecker fc Yorke 



20071 ). Since high-mass 



physical processes involved are still elusive (e.g., 
stars generally form in clusters, addressing the formation of the host cluster is a more 
fundamental process in the study of high-mass star formation, and one which would beneflt 
from more attention. 

External effect s, such as multiple sup ernova explosions that form supershells (e.g.. 



Dawson et al. 



2OIII ) and HII regions (e.g.. 



Yamaguchi et al. 



200l[ ). may play an important 



role in the formation of clusters, because compressive triggering increases the external 
pressure, promoting massive cluster formation. Similarly, in the Antennae galaxies, a 
galaxy-galaxy collision played a role in triggering clust er formation on kp c scales, and 



2000). The collision 



several young globular clusters may have been formed (I Wilson et al. 
must have created strong shocks with a relative speed of ~100 km s~^ or more, which 
presumably lead to the rapid formation of the natal massive and compact pre-cluster 
clouds. The distance of these galaxies ~20 Mpc, however, does not allow one to resolve 
the distribution of molecular gas toward individual clusters, and the detailed mechanism of 
cluster formation here remains as an open issue. 



More recently. 



Furukawa et al. 



(I2OO9I ) discovered two giant molecular clouds (GMCs) 
toward the young cluster Westerlund 2 that are likely to have recently collided with each 
other. Westerlund 2 is a remarkable Galactic super-star cluster which includes more than 
10 high-mass stars with total stellar mass of ~4500 Mq in a small volume only ~1 pc in 
radius, and which has an age of a few Myrs and is located at a distance of 3-8 kpc in 
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the G alaxy (lAscenso et al.l 120071 : iFurukawa et al.ll2009l : iRauw et al.l 120071 ). iFurukawa et al 



( 120091 ) presented morphological evidence for the physical association of the clouds with 
the cluster and its ionizing HII region RCW49 and argued that the cluster may have been 
formed via a collision between the t wo GMCs as suggest ed by their relative positions, 



masses and velocities. Subsequently, 



Ohama et al 



( I2OIOI ) derived significant temperature 



enhancements in the GMCs to ward RCW49, o ffering further robust verification of their 



association with Westerlund 2. 



Higuchi et al 



(|2010[ ) have also presented a study of 14 



young star clusters associated with dense molecular clumps of several hundreds of solar 
masses, and suggested clump-clump collisions as a possible formation mechanism for some 
of these objects. However, studies of the cloud-cloud collision scenario are yet limited 
to a small number of clusters and it is important to find more collision-triggered cluster 
formation in order to better understand the formation of clusters and high-mass stars. 



The Trifid Nebu: 



dust lanes (Walborn 



a fM20, NGC6514) is a we : 



I973I: 



Chaisson fc Willson 

— n 



1-kno wn open cluster with three outstanding 



19751). with an estimat ed distance from 1.7 



kpc fiLynds. Canzian fc 0'neilill985i) to 2.7 kp c fiCambresv et al 



detected 320 stars within 10 pc of the center (lOgura fc Ishida 



study identified 85 T-Tauri stars in the nebula (IRho et al 



20111) . An optical study 



I975I ). and a near- infrared 



20011 ). The total stellar mass 



of M20 is estimated to be less than 500 Mq, several times less massive than Westerlund 
2. M 20 also harbors one of the youngest known clusters with an estimated age of ~0.3 



Myr (jCernicharo et al 



19981). Ionization of the HII region, which extends by ~3-5 pc in 
diameter, is dominated by an 07.5 star HD164492, and some of the low-mass stars in the 
cluster appear to have bee n formed togeth er with this O star as suggested by their spatial 



proximity of a few pc (e.g., 



Rho et al. 



20081 ) . The mass- loss rate of the 0- star and expansiori 



veloc ity of the stellar wind are estimated to be M ~ 2 x 



19891) and ~ 1582 km s"^ f lPrinja. Barlow fc Howarth 



QzIMq yr"^ flHowarth fc Prinia 



19901), respectively, and the total 



dynamical luminosity of the wind is calculated to be 1/2MV^ ~ 1.6 x 10^^ erg s ^. In 
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addition, M20 is an on-going star forming region. Observations at multiple wavelengths 
revealed evidence for second generation star formation triggered b y the first generation 



stars and for relat ed activities including optical jets (e.g., HH 399, 



Hester et al. 



)J) mi 



20041 ) mid- and far- infrared protoste 



and X-ray young s t ellar objects ( 



stars (iHerbig 



1957 



Rho et al. 



2001 



lar objects ( 



2004 



Yusef-Zadeh. Biretta fc GebalL 



Cernicharo et al. 



Rho et al. 



jefioch et al. 



1998 



20061 ). infrared 



20011), Ha emission 



20051), etc. M20 therefore offers a good 



opportunity for a comprehensive study of cluster formation; from the first generation 0-star 
and many low mass stars to the subsequent second generation young stellar objects. 

Detecting and characterizing any residual natal molecular gas is crucial for 
understanding the form ation of clusters. Detai led studies of mole cular line emission in 



M20 were presented by 



fl2000h and 



Cernicharo et al. 



Lefloch et al. 



( Il998l ) and extended by 



Lefloch fc Cernicharo 



(120081 ). They covered a 20' x 20' area centered on the 07.5 star 
HD 164492 with a high angular resolution of few tens of arcseconds and detected molecular 
structures apparently associated with M20 in the velocity range 0-30 km s~^. Parts of 
the molecular features apparently trace the three dark lanes. They also identified cold 
dust cores around M20 from observations of millimeter continuum emission and discussed 



on-going star for mation, with candidates 



Spitzer telescope ( ICernicharo et al. 



1998 



"or young ste 


ar o 


Lefloch et al. 


2008 



ects (YSO) det ected by the 



Rho et al. 



2006|). However, 



it is not yet understood how these molecular clouds are organized on a scale covering 
the whole cluster, and the possiblility of triggered formation of the cluster has not been 
explored. 

We present here new CO J=2-l and J=l-0 observations of M20 with the 
NANTEN2/NANTEN 4m sub-mm/mm telescope in Chile. Our principle aim is to reveal 
degree-scale distribution of molecular gas in order to investigate the formation mechanisms 
of the first generation stars in M20, especially the central O star and the associated cluster. 
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This paper is organized as follows; Section 2 summarizes the observations and Section 3 the 
results. Discussion is given in Section 4 and a summary in Section 5. 



2. Observations 

Observations of the J=2-l transition of CO were made with the NANTEN2 4 m 
sub-millimeter telescope of Nagoya University at Atacama (4865 m above the sea level) in 
Chile in January- February and October- November 2008 for ^^CO and in April-November 
2008 for i^CO. The half-power beam width (HPBW) of the telescope was 90" at 230 GHz. 
The 4 K cooled SIS mixer receiver provided a typical system temperature of ~200 K in 
the single-side band at 220-230 GHz, including the atmosphere toward the zenith. The 
spectrometer was an acousto-optical spectrometer with 2048 channels, providing a velocity 
coverage of 392 km s^^ with a velocity resolution of 0.38 km s^^ at 230 GHz. We observed 
a large area of M20 and the supernovae remnant (SNR) W28 in ^^C0(J=2-1), while 
^■^C0(J=2-1) observations are limited in a smaller shown in Figure [H The OTF 

(= on-the-fly) mapping mode was used in the observations, and the output grid of the 
observations is 30". We smoothed the velocity resolution and spatial resolution to 0.95 km 
s~^ and 100", respectively, to achieve a better noise level. Finally, we obtained the rms 
noise fluctuations of ~0.3 K and ~0.1 K per channel in ^^CO and ^'^CO, respectively. 

Observations of the ^^00(7=1-0) and ^^CO{J=1-0) transitions were made with the 
NANTEN 4m millimeter telescope of Nagoya University at the Las Campanas Observatory 
in Chile and were carried out with a 4 K cryogenically co oled Nb superconductor-insulator- 



superconductor (SIS) mixer receiver (lOgawa et al 



1990) during the period from March 



1999 to September 2001 for ^^CO{J=1-0) and in October 2003 for ^^CO{J=1-0). The 
spectrometer was an acousto-optical spectrometer (AOS) with 2048 channels, which gave a 
frequency band-width and resolution of 250 MHz and 250 kHz, corresponding to a velocity 
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coverage of 650 km and a velocity resolution of 0.65 km at 115 GHz. The HPBW 
of the telescope was 2'. 6. Observations were made in the position-switching mode with a 
4' grid spacing. The typical system temperature was ~280 K in the single-side band, and 
typical rms noise level was 0.36 K at a velocity resolution of 0.65 km s~^. 



3. Results 

3.1. Observational results 

Figure [T] shows the ^^CO(J=1-0) integrated intensity distributions between V\sr of 
and 30 km s~^ in the region (/, 6) ~ (5°. 5- 7°. 5, — 0°.8-0°.7). M20 is denoted by a cross and 
is located at ~ 0°.6 from the supernova remnant (SNR) W28. The distribution of molecular 
clouds is very complicated, and includes an elongated feature running from south-west to 
north-east throughout the region. Dense molecular clouds associ ated with W28 are fo und 



around {l,b) ~ (6°. 4-6°. 8, — 0°.4-0°.0) in this velocity range (e.g.. 



Arikawa et al. 



19991) 



The region shown in Figure [2] corresponds to the box drawn in solid lines in Figure 
[U Figures |2^-[2t show the integrated intensity distributions of ^^C0(J=2-1) over three 
different velocity ranges. Figures [2li-[2]F show comparisons between the ^^C0(J=2-1) and 
the Spitzer 8m data, and Figures E^-M show comparisons between the ^^C0(J=2-1) and 
the IRAS 25 /im data. The large cro ss depicts the centr al star. Class O/I objects and cold 



dust cores (cloud cores TC00-TC17 



Leflochetal 



20081 ) are shown by circles and small 



crosses, respectively. These objects indicate recent star formation and pre-star formation of 
the second generation of stars in M20. 

In Figures [2^, [2li and [2^, we find a cloud with its peak position near the central star 
and well aligned with the center of the Trifid dust lanes (hereafter 2 km s~^ cloud). The 
distribution of the molecular gas shows a triangle shape and each of its corners appears to 
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trace the three dust lanes. 

The molecular gas around 8 km s~^, shown in Figures [2)3, [2fe and [2Ji, consists of 
a central cloud and three surrounding clouds to the northeast, northwest and south of 
(hereafter clouds C, NW, NE and S, respectively). In addition, we list three small clouds 
named clouds NEl, NE2 and NE3, respectively, as indicated in Figure |2b- The peak 
position of cloud C corresponds to that of the 2 km s~^ cloud. Cloud NW is associated with 
TCO without any Class O/I objects. Cloud S are located at south of cloud C with TC2, and 
Cloud NE appears to be located at the top of a molecular filament extending to the east 
from cloud C. 

The molecular gas around Visr ~ 18 km s~^ shown in Figures [2]F and [2} 
extends from the north to the south throughout the region, and we found that there 
is a elongated molecular feature overlapping with the central star on its north at 
{Ra,Dec) - (18'^02"10"-18'^02'"40", -23°16' — 22°46') (hereafter, 18 km s^^). Some of 
the other molecular co mponents in this veloci ty range are likely associated with a infrared 



dark filament in M20. 



Cambresy et al. 



(120 111 ) studied the extinction towards M20 and 



found that cold dust cores TCOO, TCO, TC3, TC4 and TC5 are located toward the 
dark filaments (Figure [3]). These cold cores show molecular counterparts in this velocity 
range. A cloud containing TC3, TC4 and several Class O/I objects are distributed at 
{Ra,Dec) ~ (18'^2'^6^ -23°6'0"), just west of the 18 km s'^ cloud (hereafter TC3 & TC4 
cloud). We also found molecular clouds associated with the Tri fid Junior, which is another 



infrared nebula at Dec 22°46' — 22°50' detected by Spitzer fiRho et al.ll2008h . The total 

integrated intensity of ^^CO(J=1-0) in this velocity range accounts for about 60 % of the 
sum of those in three velocity ranges, meaning that a large fraction of molecular mass in 
this region is concentrated to around Visr of ~18 km s~^. 

The detailed distribution of the molecular gas shows good positional coincidence 
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with optical dark features (Figure |3]) . The velocity ranges of the three panels in Figure 
12] correspond to those of Figures [2^, b and c. The current low angular resolution is not 
sufficient to resolve individual correspondences. However, the general tend may still be 
discerned in Figure |3l We find that many molecular features correspond well to the dark 
lanes and dark clouds in the optical image, indicating that they are located on the front 
side of, or partially embedded in, the nebula. 

Figure [3^ shows that the 2 km s~^ cloud is divided into two components and each of 
them trace the four major dust lanes elongated from the center to the east, the west, the 
north and the southeast. Cold dust cores TCI, TCll, TC8, TCIO and TC13 are associated 
with this cloud (Lefioch et al. 2008). Figure [3]d shows that cloud C coincides with the 
center of the nebula whereas it does not match well the dark lanes, suggesting that cloud 
C is located within or on the far side of the nebula. Cloud S seem s to correspond to t he 



pillar-like feature in the south of the nebula associated with TC2 ( iLefloch et al 



20021 ) 



indicating that it is located within or on the far side of the nebula. Cloud NW, which 
has the highest ^^CO intensity, corresponds well to a dark feature in the northwest, and 



i ntera ction s between cloud NW and the HIT region have been suggested by 



fll998h and 



Lefioch et al. 



Cernicharo et al. 



( I2OO8I ). Two small features in the north, clouds NEl and NE3, 
correspond to small dark lanes and cloud NE2 corresponds to the eastern extension of the 
2 km s"^ cloud along the major dark lane. Figure [St sh ows that the 1 8 km s ~^ cloud has a 



Lefioch et al 



(I2OO8I) indicates that 



peak component to the southwest of the central star and 
TC9 is associated with the cloud. The 18 km s~^ cloud has no counterpart dark features in 
the optical image and is likely lo cated on the far side o f M20. The total Ay in the vicinity 



of the central star is 10-15 mag (jCambresy et al. 



201 ll ). while a direct measurement towa rd 



19851), 



the central star itself gives a quite small Ay of 1.3 mag ( iLynds. Canzian &: O'neill 
indicating that a large fraction of the gas toward the central star is located on the far side 
of the nebula. We suggest that although the contours of the 2 km s~^ cloud and cloud 
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C appear to overlap with the O star, this is hkely the resuh of poor resolution in CO, 
which spreads the emission over a wider spatial area that it in fact occupies. The high 
resolution dust image in Figure [2Ji-f indicates that the dust filaments are located toward 
the west of the O star with some offset, suggesting that the amount of actual obscuring 
matter is minimal toward the O star. This minimal obscuration is consistent with a far-side 
location for the 18 km s~^ cloud. The TC3 & TC 4 cl oud is also clearly seen in Figure 



3c. 



Lefloch &: Cernichard (|2000( ) and 



Rho et al 



(120081 ) studied this cloud in detail and 



indicated that it is associated with the HII region at its boundary. 

Figure S] shows a position- velocity diagram of the molecular clouds. The integration 
range in Declination by the dotted lines is shown in Figures |2^-[2l:;. Three velocity 
components are clearly seen. We find that the 2 km s~^ cloud has a positive velocity 
gradient to the east, while the 8 km s~^ clouds have a negative velocity gradient to the east, 
and that these two velocity components seem to be connected with each other at Visr ~6 
km s~^ around the position of the O star. On the other hand, the 18 km s~^ clouds are 
extensively distributed and are detached from the 2 km s~^ and 8 km s^^ clouds in velocity 
space. 

We define seven clouds, the 2 km s"^ cloud, clouds C, NW, S, NE, the 18 km s~^ 
cloud and TC3 & TC4 cloud, in ^^C0(J=2-1) by plotting contours at half of the peak 
integrated intensity shown in Figure |2^-c, and if other peaks are found inside the contour, 
we distinguish these by plotting a line at the valley of the intensity distributions. (We 
here use a threshold of two-thirds of the peak intensity for the 18 km s"^ cloud and TC3 
& TC4 cloud because the contour at half peak intensity contains too many other peaks.) 
Figure |5] shows the observed spectra at the peak positions of the seven clouds. We here 
assume two heliocentric-distances, 1.7 kpc and 2.7 kpc for estimating sizes and masses of 
the clouds. Details of the identified clouds are shown in Table [TJ Here we can not measure 
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velocity information for the TC3 & TC4 cloud because of self absorption in the spectrum 
at a velocity of ~21 km s^^. Velocity widths of the spectra at the peaks are ~3-6 km s~^, 
with the 18 km s~^ cloud showing a somewhat wider spectrum than the others. The radii 
of the clouds are estimated to be ~l-4 pc. The masses are estimated to be a ~l-2xlO^M0 
for the 2 km s~^ and 8 km s~^ clouds and ~3-4xlO^M0 f or the 18 k m s~^ cloud. We 



here assume an X-factor of 2.0 x 10^° cm^^ (K km s^^)^^ ( .Strong et al.l 119881 ) to convert 
^^CO(J=1-0) integrated intensity into H2 column density in estimating these masses. As 
a check against previous work, we may estimate the entire H2 mass of this region within 
the full velocity range ( — 1.9-26.3 km s~^). This results in a derive d mass of 8.5x1 -^ M^ri 



Cambresy et al. 



(1201 ih 



and 2.2x10^ M0 for distances of 1.7 kpc and 2.7 kpc, respectively, 
estimated the mass of a somewhat larger area in this region to be 5.8x10^ M© by assuming 
a distance of 2.7 kpc, which is comparable to our estimate. 

Figure [6] shows the spatial distribution of the ratio between the ^^CO(J=1-0) and 
^^C0(J=2-1) transitions, which represents the degree of the rotational excitation of the 
molecular gas and reflects temperature and density. We find that the 2 km s~^ cloud and 8 
km s~^ clouds C, NW and S show significantly higher ratios of > 0.8, while the cloud NE 
shows somewhat smaller value of 0.6 at its maximum. Clouds at V\st ~ 18 km s^^ including 
the 18 km s~^ cloud and TC3 & TC4 cloud show typically smaller ratios of ~0.5 and up to 
~0.8. These results are also seen in Figure |5l The molecular gas at the southern edge of 
Figure |2]d is associated with the nearby SNR W28 and shows exceptionally higher ratios of 
~1.0. 



3.2. Physical parameters of the molecular clouds 



We estimate t he temperature and 



LVG analysis (e.g. 



Goldreich &: Kwan 



dens ity of the M20 molecular clouds with an 
1974 ). We took i2C0(J=2-l), ^'^CO{J=1-0) and 
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^^C0(J=2-1) line ratios at 6 positions, as shown in Figure [5] and Tabled] We did not use 
^^CO(J=1-0) because this transition is optically thick and traces mainly the lower density 
envelopes of the molecular clouds, and we did not estimate parameters for the TC3 & TC4 
cloud because it shows self absorption in its ^^CO spectra, which prevent us from measuring 
accurate line intensities. The LVG approximation assumes a cloud with uniform ve l ocity 



gradient, which may not alw ays be correc t in HII regions. However, ILefung fc Liszt 



and subsequent studies (e.g.. 



White 



fll976h 



19771 ) investigated radiative transfer calculations with 



a microturbulent cloud associated with HII regions and concluded that there are not any 
remarkable differences between the LVG and microturbulence models. We therefore adopt 
the LVG approximation in the following. 



We here assume abundance ratios of [^^C01/ i^^C01=77 (IWilson &: Rood 



1994h and X(CO) = [^^CO 



Leung. Herbst fc Huebner 



[Ho 



10- 



le.g. 



Frerking. Langer fc Wilson 



1982 



1984J ). respectively. We estimate velocity gradients of 1.6 
km s"^ pc~^ and 1.0 km s"^ pc~^, and therefore X (CO) / {dv / dr) of 6.3 x 10"^ and 1.0 x 10""^ 
(km s^^ pc^^)^^, for distances of 1.7 kpc and 2.7 kpc, respectively. We derived dv/dr by 
taking the average ratio between the cloud size and velocity width for the six clouds (Table 
[1]), and then took an average of them. Figures [7] and [8] show the results of the calculations, 
and details of the results are summarized in Table [21 The density and temperature ranges 
covered by the present analysis are 10^-5 x 10^ cm~^ and 10-100 K, respectively. In 
both figures, the 2 km s~^ cloud and clouds with higher line ratios of > 0.8, show kinetic 
temperatures of ~30-50 K and densities of ~ 10'^ cm^, where slightly higher temperatures 
are found in the D=2.7 kpc case. Temperatures of ~30-50 K are higher than the typical 
temperatures of dark clouds without a heat source, which are typically ~10 K. The cloud 
NE shows lower temperature of ~10 K. 



- 13 - 



4. Discussion 



4.1. The parent cloud(s) of the cluster 



Second generation star formation is currently ongoing in M 20, as evidenced by the 



Rho et al 



20081). On the 



many young stellar objects (YSOs) in the clouds (see review by 
other hand, the stars of the first generation were likely formed together with the central O 
star. An age of 0.3 My rs is estimated for the HII region by considering its spatial extent 



(ICernicharo et al. 



1998h and the 



1 Myrs (see review by IRho et al. 



age o f the first generation stars is thought to be less than 



20081 ). 



The distribution of the IRAS emission shown in Figures |2^-[2} is centered on the O 
star, and the bou ndary at the half peak intensity level shows an extension similar to that 



of the HII region (ICernicharo et a. 



following equations ( 



Dale fc Helou 



2OO2I); 



1998r) . The total infrared luminosity Li^. is given by the 



= 2.4031/25^^25 - 0.2454i/6oi^60 + 1.6381z/iooi^ioo [Lq] 



(1) 
(2) 



Here fx is the flux in each of the 3 IRAS bands and D is the distance to the object. Lir 
for M20 is hence estimate d to be ^ 2.4 x 10^ Lq and 6.1 x 10^ Lq for distances of 1.7 kpc 



an d 2.7 kpc, resp e ctive 
by 



■y- 



Draine&Li 



(I2OO7I ) pointed out that the above equations given 



Dale fc Heloul (120021 ) underestimate Lij. by up to 30% within a radiation intensity U 
range of 10-100, and applying this correction increases our estimates of to 3.4 x 10^ 
Lq and 8.6 x 10^ Lc?) for the two distances. The bolome tric luminosity of the O star is 



1.3-4.2 X 10^ Lq (iConti fc Alschuler 



1971 



Walborn 



19731 ). corresponding to 40-120 % and 



20-50 % of the above estimates for two distances, respectively. These figures indicate that 
the dust luminosity is dominated by the energy of the O star, which accounts for most of 
the radiative energy in M20. We therefore argue that the first generation stars, including 
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the O star, dominate the total energy release in the region and probably also make up the 
majority of the stellar mass contained within it. On the other hand, the stars formed in the 
other outer clouds NW, S and NE are much less dominant in luminosity and perhaps also 
in terms of the masses of individual stars. 

We here discuss the heating mechanism of the warm gas in the 2 km s~^ and clouds 
C, S and NW. The warm gas is distributed only within the extent of the strong infrared 
emission, except in the direction of the SNR (Figures [2^-[2}). It therefore seems reasonable 
to consider that the warm gas is heated predominantly by the O star, since the warmest 
regions of the 2 km s~^ and 8 km s~^ clouds are located within a few pc of this central 
source. The energy of the stellar wind from the central O star is estimated as 1.6 x 10^^ 
erg s~^ (see section 1). The cooling rate of a cloud with a dia nieter of 3 pc, density of 10 '^ 



19781 ). 



cm~'^ and temperature of 40 K is estimated as ~ 10^^ erg s~^ (iGoldsmith &: Langerl 
much lower than the energy of the stellar wind. The high temperatures of the 2 km s~^ and 
8 km s~^ clouds can therefore be explained energetically in terms of the shock heating by 
the stellar winds of the central star, although the present low angular resolution is not high 
enough to test this possibility further by deriving a detailed temperature distribution. In 
either case, the physical association of the four clouds with the O star is virtually certain. 

It is therefore reasonable that both the 2 km s~^ cloud and cloud C are candidates for 
the parental cloud of the O star. The two clouds share peak positions, suggesting that they 
are located nearly at the same position, and the velocity difference of these two clouds is 
~7.5 km s~^ (Tabled]), meaning that they can only move by ~2 pc within the 0.3 Myr age 
of M20, which is consistent with the interpretation that these two clouds is located at the 
same position. 

On the other hand, the molecular clouds at Visr ^ 18 km s~^ do not seem to be directly 
associated with the O star, because these clouds have temperatures of 10 K, lower than 
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those of the 2 km and 8 km clouds. The clouds around Visr ~ 18 km s~^ therefore 
appear to be located near M20 but are neither as close or as heated by M20 as the 2 km 
s~^ and 8 km s~^ clouds. Only TC3 & TC4 cloud has a possibility to be located arou nd 
M20 as close as cloud NE, NW and S, because studies by 



Rho et al. 



Lefloch fc Cernichard ( 120001 ) and 



(l2008l ) indicate that it is located at the boundary of the HII region and shocks 
may occur there. Self absorption in the ^^CO spectra possibly make it difficult to know an 
actual excitation condition of the cloud, and further study in optically-thin lines is needed 
for better understanding. 



4.2. Possible triggering of the cluster formation 

The masses of the 2 km s~^ cloud and cloud C, assuming distance of between 1.7 and 
2.7 kpc, are estimated to ~ 0.8-2.0 x 10^ Mq and ~ 0.5-1.2 x 10^ Mq, respectively, giving 
a total mass of about ~ 1.3-3.2 x 10^ Mq. We shall first discuss the gravitational condition 
of these two clouds. These clouds, both peaking toward the O star, have radii of ~l-2 pc 
and a velocity difference of ~7.5 km s^^ (Table [1]). Based on these figures we estimate 
that a mass of ~1.3-2.6 x 10^ Mq is required to gravitationally bind the two clouds for 
an assumed separation 1-2 pc. This mass is apparently larger than the observed cloud 
and stellar masses, suggesting that these two clouds are not gravitationally bound to M20. 
Mechanical expansion therefore may be reasonable explanation for the velocity difference. 
Since M20 is a very young object with an age of 0.3 Myr, it does not yet appear to have 
hosted any supernova explosions, and stellar winds are the only available mechanism of 
driving any mechanical expansion that may be occurring around the central cluster. The 
total mechanical luminosity available from stellar winds of the O star over the age of the 
system is estimated to be ~ 1.5 x 10^^ erg. Assuming the velocity separation of ~7.5 km s~^ 
as the expanding velocity of the clouds, we obtain the total kinetic energy of the five clouds. 
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2 km cloud, Cloud C, NW, S and NE, of ^2.0-4.7x10"^^ erg, which corresponds to 
almost 15-30 % of the energy given from the stellar winds. Numerical calculations suggest 
that le ss than a few percent of the initial wind energy is converted to neutral gas kinetic 



energy ( lArthur 



20071 ). and the small sohd angle less than 4 vr covered by the clouds reduces 
the energy available to them. We therefore conclude that while mechanical expansion 
cannot be ruled out as the cause of the velocity separation, the energetics are uncomfortably 
tight. Moreover, the spatio-velocity structure of the two components, as shown in figure HI 
shows the greatest departures from the centroid velocity at around 18'*2"*00* - the furthest 
projected distance from the central star. This does not fit the classical pattern of a shell-like 
expansion around a central source, and suggests that at least some component of the cloud 
velocity separation is unrelated to their interaction with the central star 

We here suggest that much of the observed motion of the clouds does not arise as 
a direct result of the stellar cluster but is instead systemic — i.e. present from before 
the cluster's birth. Yet the 2 km s~^ cloud and cloud C are clearly associated with M20 
and in close spatial proximity to one another. We suggest that a scenario in which a 
cloud-cloud collision triggered the formation of the central stars is highly consistent with 
these observational characteristics. In this scenario the 2 km s~^ cloud and cloud C collided 
each other ~1 Myr ago at a relative velocity of ~7.5 km s~^ or more. The rapid collision 
between the two clouds strongly compressed the molecular gas and triggered the formation 
of the central star and surrounding first generation stars. In this model, the 2 km s~^ 
cloud and cloud C must be moving in the opposite directions; the former is moving toward 
us and the latter is moving away from us. As shown in Section 3.1, the 2 km s^^ cloud 
is apparently located at the front side of M20, whereas cloud C is not in the near side, 
suggesting that it is located either on the far side or within M20. This relative configuration 
is in fact consistent with the cloud-cloud collision scenario, because we must be witnessing 
a moment after the collision occurred Myr ago. We would expect a reversed cloud location 
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prior to the collision. 



T 

on by 



lis scenario is simila r to that discussed by 



Furukawa et al. 



Ohama et al. 



( 120091 ) and later expanded 



(120101 ). These authors found two GMCs closely associated with the 
RCW 49 HII region and its exciting cluster, the super star cluster Westerlund 2. The two 
clouds have a velocity difference of 15 km s~^, are not gravitationally bound to one another, 
and have associated kinetic energies that are too large to be explained by expansion driven 
by the central cluster. These characteristics are quite similar to the 2 km s~^ and 8 km 
clouds of M20, and we therefore suggest that M20 may be interpreted as a miniature of 
Westerlund 2. 

The velocity of the M20 re gion as a who le (0-30 km s"^) corresponds to the velocities of 



the Perseus and Scutum Arms (IVallee 



in the region (see Figure [T] and 



20081). and the re are actually many molecular clouds 



Takeuchi et al 



2010l ). Such a high density of molecular 



clouds implies a high cloud-cloud collision probabilit y, providing a s i milar environment to 



Dobashi et al. 



( I2OOII ) investigated the 



that in the region of Westerlund 2. We also note that 
relationship between masses of star forming molecular clouds and the infrared luminosity 
from these clouds, finding that molecular clouds with masses of ~ 10^-10^ Mq typically 
have infrared luminosity of up to ~ 10^-10^ Lq — almost comparative to or much lower 
than the infrared luminosity of M20 (~ 2.4-6.1 x 10^ -^0)- This fact suggests that the 
first generation star formation in M20 is remarkably efficient, providing further tentative 
support for the triggered star formati on mechanism suggested here. 



5. Summary 



We summarize the present works as follows; 

1) A large-scale study of the molecular clouds toward the Trifid nebula, M20, has 
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been made in the J— 2-1 and J— 1-0 transitions of ^^CO and ^^CO. These observations 
reveal two molecular components both peaked toward the center of M20 with LVG-derived 
temperatures of 30-50 K — significantly higher than the 10 K of their surroundings. 

2) The close association of the clouds with the central cluster and surrounding HIT 
region strongly suggests that they are the parent clouds of the first generation stars in M20. 
The mass of each cloud is estimated to be ~ 10^ Mq and their velocity separation is ~ 7.5 
km s~^ over a length of ~l-2 pc. The total stellar and molecular mass is too small by an 
order of magnitude to gravitationally bind the system. 

3) Based on the dynamics and energetics of the system, we suggest that the formation 
of the first generation stars, including the central ionizing 07.5 star, was triggered by the 
collision between the two clouds on a short time scale of ~1 Myr. This is a second case 
alongside that of the super star cluster Wester lund 2, for which a similar collision-triggered 
formation mechanism has been suggested. 

NANTEN2 is an international collaboration of ten universities, Nagoya University, 
Osaka Prefecture University, University of Cologne, University of Bonn, Seoul National 
University, University of Chile, University of New South Wales, Macquarie University, 
University of Sydney and Zurich Technical University. NANTEN was operated based on a 
mutual agreement between Nagoya University and the Carnegie Institution of Washington. 
We also acknowledge that the operation of NANTEN can be realized by contributions 
from many Japanese public donators and companies. The work is financially supported 
by a grant-in-aid for Scientific Research (KAKENHI, no. 15071203, no. 21253003, and 
no. 20244014) from MEXT (the Ministry of Education, Culture, Sports, Science and 
Technology of Japan) and JSPS (Japan Society for the Promotion of Science) as well as 
JSPS core-to- core program (no. 17004). This research was supported by the grant-in-aid 
for Nagoya University Global COE Program, Quest for Pun damental Principles in the 



- 19 - 



Universe: from Particles to the Solar System and the Cosmos, from MEXT. Also, the work 
makes use of archive data acquired with IRAS and Spitzer data sets gained with Infrared 
Processing and Analysis Center (IPAC). 



-20- 
REFERENCES 

Arthur, S. J., 2007, in lAU Symp. 250, Massive Stars as Cosmic Engines, ed. F. Bresolin, 
P. A. Crowther, & J. Puis (Cambridge: Cambridge Univ. Press), 355 

Arikawa, Y., Tatematsu, K., Sekimoto, Y., & Takahashi, T., 1999, PASJ, 51, 7 

Ascenso, J., Alves, J., Beletsky, Y., & Lago, M. T. V. T., 2007, A&A, 466, 137 

Cambresy, L., Rho, J., Marshall, D. J., & Reach, W. T., A&A, 527, 141 

Cernicharo, J. et al., 1998, Science, 282, 462 

Chaisson, E. J. & Willson, R. F., 1975, ApJ, 199, 647 

Conti, P. S., & Alschuler, W. R., 1971, ApJ, 170, 325 

Dawson, J. R., McClure-Crifhths, N. M., Kawamura, A., Mizuno, N., Onishi, T., Mizuno, A., 
& Fukui, Y., 2011, ApJ, 728, 127 

Dale, A. D., & Helou, C, 2002, ApJ, 576, 159 

Dobashi, K., Yonekura, Y., Matsumoto, T., Momose, M., Sato, F., Bernard, J., & 
Ogawa, H., 2001, PASJ, 53, 85 

Draine, B. T., & Li, A., 2007, ApJ, 657, 810 

Frerking, M. A., Langer, W. D., & Wilson, R. W., 1982, ApJ, 262, 590 

Furukawa, N., Dawson, J. R., Ohama, A., Kawamura, A., Mizuno, N., Onishi, T., & 
Fukui, Y., 2009, ApJ, 696L, 115 

Goldreich, P., & Kwan, J., 1974, ApJ, 189, 441 

Glodsmith, P. F., & Langer, W. D., 1978, ApJ, 222, 881 



Herbig, G. H., 1957, ApJ, 125, 654 

Hester, J. J., Desch, S. J., Healy, K. R., & Leshin, L. A., 1990, Science, 304, 1116 
Higuchi, A. E., Kurono, Y., Saito, M., & Kawabe, R., 2010, ApJ, 719, 1813 
Howarth, 1. D. & Prinja, R. K., 1989, 69, 527 
Lefloch, B. & Cernicharo, J., 2000, ApJ, 545, 340 

Lefloch, B., Cernicharo, J., Cesarsky, D., Demyk, K., & Rodriguez, L. F., 2001, A&A, 368, 
13 

Lefloch, B., Cernicharo, J., Rodriguez, L. F., Miville-Deschenes, M. A., Cesarsky, D., & 
Heras, A., 2002, ApJ, 581, 335 

Lefloch, B., Cernicharo, J., & Pardo, J. R., 2008, A&A, 489, 157 

Leung, C. M., Herbst, E., & Huebner, W. F., 1984, ApJS, 56, 231 

Leung, C. -M., & Liszt, H. S., 1976, ApJ, 208, 732 

Lynds, B. T., Canzian, B. J. & O'neil, E. J., Jr., 1985, ApJ, 288, 164 

Ohama, A. et al., 2010, ApJ, 709, 975 

Ogawa, H., Mizuno, A., Hoko, H., Ishikawa, H., & Fukui, Y., 1990, Int. J. Infrared and 
Millimeter Waves, 11, 717 

Ogura, K., & Ishida, K., 1975, PASJ, 27, 119 

Prinja, R. K., Barlow, M. J., Howarth, I. D., 1990, ApJ, 361, 607 

Rauw, G. et al, 2007, ARA&A, 463, 981 

Rho, J., Corcoran, M. F., Chu, Y., & Reach, W. T., 2001, ApJ, 562, 446 



- 22 - 

Rho, J., Ramirez, S. V., Corcoran, M.F., Hamaguchi, K., & Lefloch, B., 2004, ApJ, 607, 904 

Rho, J., Reach, W. T., Lefloch, B., & Fazio, G. G., 2006, ApJ, 643, 965 

Rho, J., Lefloch, B., Reach, W. T., & Cernicharo, J., 2008, in Handbook of Star Forming 
Regions, vol II, ASP Conf. Ser., 5, 509 

Strong, A. W. et al, 1988, A&A,207,1 

Takeuchi, T. et al., 2010, PASJ, 62, 557 

Vallee, J. P., 2008, AJ, 135, 1301 

Walborn, N. R., 1973, AJ, 78, 1067 

Wilson, T. L., & Rood, R., 1994, ARA&A, 32, 191 

Wilson, C. D., Scoville, N., Madden, S. C., Charmandaris, V., 2000, ApJ, 542, 120 
White, R. E., 1977, ApJ,211,744 

Yamaguchi, R., Mizuno, N., Onishi, T., Mizuno, A., & Fukui, Y., 2001, PASJ, 53, 959 
Yusef-Zadeh, F., Biretta, J., & Geballe, T. R. 2005, AJ, 130, 1171 
Zinnecker, H., & Yorke, H. W., 2007, ARA&A, 45, 481 



This manuscript was prepared with the A AS macros v5.2. 



- 23 - 




Galactic Longitude [degree) 
Conrour leveJs; mliTh, IB K kmfc, Interva) leKknVa 



Fig. 1. — Integrated intensity distribution of ^^C0(J=2-1) emission towards (/, b) ~ (5.5°- 
7.5°, — 0.8°-0.7°). Contours are drawn every 16 K km s^^ from 18 K km s^^. The cross 
indicates M20, and solid lines indicates the region shown in Figure [21 The dotted circle 
indicates the boundary of radio emission from W28, and dashed lines indicate the region 
observed in ^^C0{J=2-1). 
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Fig. 2. — (a-c) Integrated intensity distribution of ^^C0(J=2-1) emission. Contours are 
drawn every 4 K km s~^ from 6 K km s~ ^. Filled circles depict Class O/I objects identified 



by Spitzer observations (IRho et al. 



20081 ). and small crosses d epict cold dust cores (TCOO- 



TC17) identified by radio observations by 



Lefloch et al. 



( I2OO8I ). The large cross depicts the 



central O star, (d-f) Contour rn aps of the ^^C0(J=2-1) emission superposed on the Spitzer 



8 /im image of 



Rho et al. 



(120061 ). Contours are drawn every 8 K km s~^ from 14 K km s~^. 
(g-i) Contour maps of the ^^C0(J=2-1) emission superposed on IRAS 25 /im data. White 
contours mark the two-third peak intensity of the IRAS emission equal to 60 MJy str^^. 
Solid contours show ^^C0(J=2-1) emission and are drawn every 8 K km s"^ from 14 K km 
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Fig. 3. — C0( J=2-l) distribution of the molecular clouds associated with M20, superposed 
on an optical image (Credit: NOAO). Contours for all clouds are drawn every 4 K km s~^ 
from 14 K km s~^. Small crosses depict cold dust cores, with ID numbers labeled only in 
Figure [3^.. The large crosses depict the central star. 
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Fig. 4. — Position-velocity diagram of ^^C0(J=2-1) emission. The integration range is 
shown in Figure l(a-c) by dashed hnes. The dashed hne in the figure shows the position of 
the central O star of M20. 
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Fig. 5. — CO spectra at the peak positions of the clouds hsted in tabled] ^^CO(J=1-0), 
12C0(J=2-1), 13CO(J=1-0) and ^^C0{J=2-1) are plotted in black, red, blue and green, 
respectively. All spectra were smoothed to be a beam size of 2. '6. Dotted lines indicated the 
velocity ranges used for the LVG calculations shown in figures [7] and [3 
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Fig. 6. — Spatial distribution of ^^C0(J=2-1) emission (contours) and ratio between 
^^C0(J=2-1) and ^^CO(J=1-0). The data used for the calculation of ratios was spatially 
smoothed with a gaussian function to an effective beam size of 2'. 6. C ontours are plott ed at 



every 5 K km s ^ from 11 K km s ^. Circles depict Class O/I objects (iRho et al. 
the cross depicts the central O star of M20. 



20081 ) and 
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Fig. 7. — LVG results for X(CO)/dvdr = 6.3 x 10^^ (km s^^ pc^^)^^, assuming a distance 
of 1.7 kpc, are shown in the density-temperature plane. Crosses denote the points of min- 
imum x^i ci-nd filled areas surrounding the crossses indicate of 3.84, which corresponds 
to the 95 % confidence level of the distributions with 1 degree of freedom. Solid, dot- 
ted and dashed lines show ^^CO{J=l-Q) /^^C0{J=2-l), ^^C0{J=2-l) /^^C0{J=2-l) and 
13CO(J=2-1)/13CO(J=1-0) intensity ratios. 
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Fig. 8. — Same as figure [7] for X(CO)/dvdr = 1.0 x 10 (km s ^ pc ^) ^, assuming a 



distance of 2.7 kpc. 
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Table 1: Physical properties of the clouds. 





Ra(J2000) 


Dec(J2000) 




AV 




M(H2)t 


Name 














s 


(° ' ") 


(km s~^) 


(km s~^) 


(pc) 


(xlO^ Mq) 


2 km s~^ cloud 


18:02:19.50 


-23:01:40.4 


1.4 


4.0 


1.2/2.0 


0.9/2.0 


cloud C 


18:02:18.36 


-23:01:45.3 


8.9 


3.6 


0.9/1.5 


0.6/1.2 


cloud NW 


18:02:04.31 


-22:56:48.8 


10.1 


3.3 


1.2/1.9 


0.9/2.3 


cloud S 


18:02:25.22 


-23:05:42.9 


7.5 


3.7 


0.8/1.2 


0.6/1.0 


cloud NE 


18:02:51.35 


-22:58.31.7 


9.9 


2.6 


1.2/2.0 


0.7/1.8 


18 km cloud 


18:02:32.31 


-23:03:25.1 


18.2 


6.4 


2.3/3.7 


5.6/14.1 


TC3 & TC4 cloud 


18:02:16.30 


-23:05:25.1 


_t 


_t 


1.0/1.6 


0.8/2.2 



'^Sizes and masses of the clouds were calculated with the both distances of 1.7kpc/2.7kpc. 
^VisT and AV of TC3 & TC4 cloud can not be measured due to self absorption in the spectrum. 

Note. — Column (1): Name of cloud. (2-3): Peak position of cloud. (4): Velocity of cloud. (5): Velocity 
dispersion (FWHM). (6): Radius of cloud. (7): Cloud mass derived using an X-factor of 2.0 x 10^°cm~^ (K 
km s-i)-^ 
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Table 2: LVG results. 



Name 






1.7 kpc 






D = 


2.7 kpc 






















n(R2) ( 


cm~^) 


71c (K) 


min. 


nm ( 


cm~^) 


rk(K) 


min. 


2 km s~^ cloud 


1 9+0.3 
-■-•^-0.3 


X 10=^ 


oq+28 


0.002 




X 10^ 


CO +40 
^-^-32 


0.006 


cloud C 


^- ' -0.2 


X 10^ 




0.008 




U-J-O.l 


X 10^ 


> 40_29 


0.0007 


cloud NW 


^■0-0.5 


X 10^ 


301}^ 


0.02 


-■-■^-0.4 


X 10^ 


381}? 


0.01 


cloud S 


1 4+0-5 
-■-■^-0.4 


X 10^ 




0.004 


1 -1+0.5 
^■^-0.3 


X 10^ 


481^° 


0.004 


cloud NE 


i.olSi 


X 10^ 


< 11+9 


0.001 


0.8l°i 


X 10^ 


< 12+11 


0.006 


18 km s~^ 




X 10^ 


< 10+1 


0.15 


0.71°:} 


X 10^ 


< 10+2 


0.08 



Note. — Column (1): Name of cloud. (2 4): LVG result for XiCO){dv/dr) = 6.3 x IQ-'"^ (km s^i pc^i) 
at a distance D of 1.7 kpc. (5 7): LVG result for X{GO){dv/dr) = 1.0 x 10"'' (km s^^ pc^^) at a distance 
D of 2.7 kpc. (2, 5): Number density of H2. (3, 6): Kinetic temperature. (4, 7) minimum x^- 



